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o 

(N 
X 



Gao-Chan Yong x 'B and Bao-An Li 2 <Q 

1 Institute of Modern Physics, Chinese Academy of Sciences, Lanzhou 730000, China 
2 Department of Physics and Astronomy, Texas A&M University-Commerce, Commerce, TX 75429, USA 
and Department of Applied Physics, Xi 'an Jiaotong University, Xi'an 710049, China 

Using a relativistic transport model ART1.0, we explore effects of nuclear symmetry energy on 
rj meson production and its rare decay to the dark U-boson in heavy-ion reactions from 0.2 to 10 
GeV/nucleon available at several current and future facilities. The yield of n mesons at sub-threshold 
energies is found to be very sensitive to the density dependence of nuclear symmetry energy. Above 
a beam energy of about 5 GeV/nucleon in Au+Au reactions, the sensitivity to symmetry energy 
disappears. Using the branching ratio of the rare n decay (n — > jU) available in the literature, we 
estimate the maximum cross section for the U-boson production in the energy range considered, 
providing a useful reference for future U-boson search using heavy-ion reactions. 

PACS numbers: 25.75.-q, 21.65. Mn, 21.65.Ef, 95.35.+d 
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I. INTRODUCTION 

The physics motivations of this work are twofold. The 
first purpose is to look for new and possibly more sensi- 
tive probes of the high-density behavior of nuclear sym- 
metry energy in high energy heavy-ion collisions. The 
other purpose is to explore the possibility of producing 
the neutral vector U-boson (dark photon) introduced in 
the super-symmetric extension of the Standard Model 
(SM) 1-5] in relativistic heavy- ion collisions. The two 
purposes naturally come together in studying properties 
of neutron stars where both the possible existence of the 
U-boson and the high-density behavior of nuclear sym- 
metry energy affect significantly the Equation of State 
(EOS) of dense neutron-rich nucleonic matter. In fact, it 
has been shown already that massive neutron stars can 
be supported even with a rather soft EOS for symmetric 
nuclear matter and/or a super-soft symmetry energy if 
the additional interaction between two nucleons due to 
the exchange of a U-boson is considered @, 0] • Moreover, 
both the symmetry energy and the possible U-boson af- 
fect the core-crust transition density/pressure and the 
moment of inertia of neutron stars [8J . 



A. Potentials of using n meson production in 
heavy-ion collisions as a probe of the high-density 
behavior of nuclear symmetry energy 

Nuclear symmetry energy E sym (p) encodes the energy 
cost of converting all protons into neutrons in symmet- 
ric nuclear matter. Because the density dependence of 
E S ym{p) is very important for understanding proper- 
ties of rare isotopes and neutron stars as well as the 
dynamics of heavy-ion reactions, supernova explosions 



and gravitation wave emissions of spiraling neutron star 
binaries, much efforts have been devoted to extract- 
ing the E sym (p) using data from both terrestrial labo- 
ratory experiments and astrophysical observations, see 
e.g., refs. [9l-l22i|. While significant progress has been 
made recently in constraining the -E S ym(po) and the slope 



L(Po) 
refs. |2 



[3p(dE sym /dp] at normal density po, see, e.g., 
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49j, our knowledge about the high-density be- 
havior of the E sym (p) is still very poor. In fact, findings 
about the high-density E sym (p) from comparing heavy- 
ion reaction data with various transport model calcula- 
tions are still inconsistent, see, e.g., refs. (50rl53| . This is 
at least partially because of both the uncertainties in the 
physics inputs and the different numerical techniques in 
initializing and modeling the transport process of collid- 
ing nuclei, see, e.g., refs. [53T - l58j for discussions on some 
of the contributing factors. Thus, new and possibly more 
sensitive probes of the high-density behavior of nuclear 
symmetry energy are always useful. Since the symme- 
try (isovector) potential is normally much smaller than 
the isoscalar potential, in order to observe clearly effects 
of the symmetry energy/potential in heavy- ion reactions 
one needs ideally to use slowly moving particles to be 
acted on by the isovector potential for a long time in a 
large region of high isospin asymmetry and density gra- 
dients. Sub-threshold mesons, such as pions and kaons, 
especially the ratios of different charge states, are good 
candidates and have been studied extensively. Transport 
model calculations have shown consistently that near 
their respective production threshold, the tt~ /ir + [59T463] ] 
and K + /K° [64 1 ratios are indeed sensitive to the high- 
density symmetry energy although comparisons with lim- 
ited data available are still inconclusive so far. Generally 
speaking, more massive mesons probe the E sym (p) at 
higher densities if effects of the final state interactions do 
not wash out the signal. The r\ meson of mass 547.853 
MeV/c 2 [HI is the most massive member in the octet of 
pseudoscalar Goldstone mesons. It has significant photon 
and dilepton decay channels providing the possibility of 
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studying the high-density symmetry energy more cleanly 
using the electromagnetic probes. In fact, 77 meson pro- 
duction in heavy-ion collisions has been studied exten- 
sively both theoretically and experimentally, especially 
by the TAPS and HADES Collaboratio ns, see, e.g., refs. 
[66M75I ] . Moreover, as pointed out earlier [Sa,[£Z|, because 
of the hidden strangeness (the ss component), 77 mesons 
experience weaker final state interactions compared to 
pions. Contrary to kaons, because the net strangeness 
content is zero in 77 mesons, they can be produced with- 
out another strange partner in the final state and thus 
require less energies. More quantitatively, the threshold 
invariant energy for 77 and kaon production in nucleon- 
nucleon collisions is (v 7 *)^ = 2.4 GeV and {yfs)^ h = 2.6 
GeV, respectively. The corresponding threshold beam 
energy in free nucleon-nucleon collisions is E^ h = 1.2 
GeV and Ef h =1.6 GeV, respectively. For a compar- 
ison, the experimental K + /ir + and 77/ 7r° ratios in the 
most central Au+Au reactions at 1 GeV/nucleon is ap- 
proximately 2 x 10~ 3 [zH and 1.4 x 10~ 2 [73], respectively. 
It is thus interesting to know if the 77 meson might be use- 
ful for exploring the E sym (p) . Intuitively, in order for 
the p-n pair to create a 77 near its production threshold, 
they have to have their momenta oriented opposite to 
each other. Thus, near threshold 77 mesons will probe 
sensitively the p-n relative momentum which is deter- 
mined mainly by the gradient of the isovector potential. 
Compared to pions and kaons, however, the elementary 
77 production cross sections in baryon-baryon and meson- 
baryon scattering still suffer from relatively larger un- 
certainties although they are gradually better known as 
more data and calculations become available [H, [78| . It 
is thus necessary to know the optimal beam energies and 
the kinematic regions where the yield of 77 mesons is most 
sensitive to the E sym (p) . 



B. Heavy- ion reactions as a possible tool for dark 
U-boson search 



While there are well established observational evi- 
dences for dark matter through its gravitational inter- 
actions, the mass and interactions of dark matter par- 
ticles are completely unknown [7j|. Thus, what is dark 
matter? It is actually at the very top of the Eleven Sci- 
ence Questions for the New Century identified by the 
Committee on the Physics of the Universe, US National 
Research Council [80( | . In fact, much efforts have been 
devoted to searching for dark matter candidates, see, 
e.g., refs. [HI, HI] for reviews. Besides the weakly in- 
teracting massive particles as the most popular candi- 
dates for dark matter, it has been proposed that light 
dark matter particles v with a mass in the range of 0.5- 
20 MeV may exist [lHa]- I n particular, the annihilation 
°f XX ~~ ^ e + e~ through the exchange of the light vec- 
tor U-boson of mass 10-100 MeV, has been used to ex- 
plain successfully the INTEGRAL satellite observation 
of the excess flux of 511 keV photons from the center re- 



gion of our galaxy [83| although alternative explanations 
exist 8Jj. Moreover, assuming the U-boson couples to 
quarks as well as electrons [79(, the extra contribution 
to the pion decay (7r° — > e + e~) width mediated by an 
off-shell U-boson can explain nicely the enhanced pion 
decay observed by the KTeV Collaboration I85| in com- 
parison with the Standard Model prediction 86]. Fur- 
thermore, as mentioned earlier, the extra interaction be- 
tween nucleons due to the U-boson exchange may affect 
properties of neutron stars [53 although it has no 
effect on properties of finite nuclei [88| . It may also cause 
the deviation from the inverse-square-law of gravity, see, 
e.g., [891 - 193 ] for recent reviews. The aforementioned stud- 
ies and findings have motivated considerable interest in 
searching for the U-boson [95| . While various upper lim- 
its on the strength and interaction range of the U-boson 
have been put forward at various length scales without 
contradicting known physics principles and existing ex- 
perimental/observational data, see, e.g., ref. [88| for the 
latest review, most of the constraints on properties of 
the U-boson are indirect. It is thus very interesting to 
see that some direct constraints on this kind of gauge 
bosons through neutral meson decays in neutrino experi- 
ments at CERN have been reported recently 96j. More- 
over, several proposals to investigate directly properties 
of the U-boson in low energy (E c . m < 10 GeV) exper- 
iments have been put forward. These include experi- 
ments usin g e + e~ colliders through the e + e~ — > jU pro- 
cess [1, 195l. l97l - [99j , meson rare decays at meson factories 
and fixed target electron- nucleus scatterings [95j. Many 
mesons can decay into the U-boson with branching ratio 
BR{X -> Y + U) « e 2 BR(X -> Y + 7) where e is the 
strength ratio of the U-SM particle to 7-SM particle cou- 
pling. This is then followed by the U decay into dileptons, 
i.e., U — > £ + £~ . Using the meson summary tables of the 
Particle Data Group, Reece and Wang estimated the U- 
producing branching ratios for several mesons 95]. The 
77 rare decay was shown to be the most promising one. 
More quantitatively, the ratio of 7/ — s- *yU, oj —> ir°U, 
<f> — > r)U and — > jU is found to be approximately 
320/225/1.3/4.4 x 10~ 3 . It is interesting to note that 
the 7 + p — > 77 + p reaction with a cross section of ap- 
proximately 70 nb will b e us ed for the U-boson search 
at the Jlab Eta Factory [lOOj ]. As we shall estimate in 
this work, the 77 production cross section in a Au+Au 
collisions from 0.2 to 10 GeV/A beam energy range is 
on the order of 60-6000 mb. While it will be very chal- 
lenging to extract information about the U-boson from 
various backgrounds, the relatively large 77 production 
cross section warrants a U-boson search using heavy-ion 
facilities, such as the CSR (Cooler Storage Ring) and the 
planned HIAF (High Current Heavy-Ion Accelerator Fa- 
cility) in Lanzhou, the FAIR (Facility for Antiproton and 
Ion Research) in Darmstadt, NICA (Nuclotron based Ion 
Collider Facility) at JINR Dubna and the BES (Beam 
Energy Scan) program at RHIC (Relativistic Heavy-Ion 
Collider). To explore the idea more quantitatively, we 
shall investigate the excitation function of 77 production 
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in Au+Au reactions in the energy range covered by these 
facilities. 



II. TRANSPORT MODEL SIMULATIONS OF 77 
MESON PRODUCTION IN HEAVY-ION 
COLLISIONS 

Our study is based on the relativistic trans- 
port model ART |lOl| . The mo del was developed 
from the BUU transport model [l02j by including 
more baryon and meson resonances and their in- 
teractions. More specifically, the following baryons 
N, A(1232), A*(1440), A*(1535), A, E, and mesons 
7r, p, oj, 77, K with their explicit isospin degrees of 
freedom are included. The model was used success- 
fully in studying many features of heavy-ion reactions 
at AGS energies up to a beam mom entum of about 15 
GeV/c, for a review, see, ref. [1031 ]. An extended ver- 
sion of ART is also used as a hadronic afterburner in 
the AMPT (A Multi-Phase Transport Model) model for 
heavy-ion collisions at RHIC and LHC energies |l04l | . In 
the ART model, rj meson production and absorption are 
modeled via r]N O A* (1535). The cross sections for 
the production and absorption of A* (1535) resonance 
in baryon-baryon c ollisi ons and its decay width can be 
found in refs. (6^ llOlj . For ease of the following dis- 
cussions, we emphasize that a(NN — » AA*(1535)) « 
2a (NN — > NNr/) and the elementary cross section for 
77 production is strongly isospin dependent, i.e., a(pn — > 
pnrj) ss 2>a{pp — > ppri) = 3a(nn nnq). The ART1.0 
has the option of using different isoscalar potentials for 
baryons. However, the isovector (symmetry) potential 
was not implemented earlier. For this exploratory study, 
we use the Skyrme-type parametrization for the isospin- 
depend ent but m omentum-independent mean field po- 
tential dM Eol 

U( P ,S,t) = A(p/p a ) + B(p/ Po y 

+ 4r£2£(p) + (18.6-F(aO) 

x (G(x) - l)(p/po) G{x) 6 2 . (1) 

In the above, EP£ n (p) = F(x)p/p + (18.6 - 
F{x))(p/ po) G ^) is the interaction part of nuclear sym- 
metry energy, 6 = (p n — Pp)/p is the isospin asymme- 
try at density p and r=l/2(— 1/2) for neutrons (pro- 
tons). Adding the kinetic part of nuclear symmetry en- 
ergy Eg^(p) = ft 2 /6m„(37T 2 p/2) 2 / 3 where m n is the nu- 
cleon mass, the symmetry energy corresponding to the 
single-particle potential in Eq. [T] is 

E sym {p) = E% n m {p) + Eff m (p) 
= n 2 /Qm n {in 2 p/2) 2 ^ 
+ F(x)p/ p + (18.6 -F(x))(p/p ) G ^ (2) 

with F(x) and G{x) given in Ref. jl05| for different values 
of the parameter x used to vary the density dependence of 



Esym [106]. Shown in Fig. [T]is the E sym (p) with x = 1,0 
and —2, respectively, from being very soft to very stiff. 
The mean-field potentials of mesons and baryon reso- 
nances in nuclear matter are still largely unknown. We 
adopt here the minimum assumption that the isoscalar 
part of the mean-field potential for baryon resonances is 
the same as that for nucleons and there is no mean-field 
for mesons as assumed in most transport models. The 
isovector potential for baryon resonances is taken as an 
average of that for neutrons and protons with an isospin- 
dependent weighting factor determined by the square of 
the Clebsch-Gordan coefficients in the A(A*) «-» n + N 
coupling [59( . Contributions of the different charge states 
of baryon resonances to the isospin asymmetry of the 
baryonic system are similarly calculated [H^] . We choose 
two sets of parameters A and B, i.e., A= -356 (-124) 
MeV, B= 303 (70.5) MeV and cr= 7/6 (2) corresponding 
to the incompressibility of nuclear matter 200 (380) 
MeV at pn = 0.168 fm~ 3 . 




FIG. 1: (Color online) Nuclear symmetry energy as a function 
of reduced density p/po (the symmetry energy at po is taken 
to be 31.6 MeV ) with three values for the parameter x. Taken 
from Ref. fw§ . 



III. EFFECTS OF NUCLEAR SYMMETRY 
ENERGY ON tj MESON PRODUCTION 

Shown in Fig. [2] is a comparison of the ART1.0 model 
calculations with the TAPS experimental data on the 
77 transverse momentum distribution at midrapidity in 
Au+Au reaction at a beam energy of 1 GeV/nucleon. 
In this and the following model calculations, the error 
bars are statistical in nature. It is seen that our cal- 
culations give a reasonable description of the data with 
x = and 1 while the stiff nuclear symmetry energy with 
x= -2 leads to a significantly larger number of 77 mesons. 
For a comparison, we notice that the FOPI ir~ /ir + data 
in Au+Au reactions at beam energies from about 0.4 
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FIG. 2: (Color online) Transverse momentum distributions of 
midrapidity r\ mesons in inclusive Au+Au reaction at a beam 
energy of 1 GeV/nucleon with different nuclear symmetry en- 
ergies. Data are taken from Ref. [70l ]. 



FIG. 3: (Color online) Multiplicity of inclusive r\ production 
as a function of incident beam energy in Au+Au reactions 
with the two different values for the symmetry energy param- 
eter x. 



to 1 GeV/nucleon favor the symmetry energy with x=l 
50]. Because 77 is sensitive to the number of p-n colli- 
sions while the tt~ /ir + ratio is determined by the ratio 
of n-ii and p-p colliding pairs, r\ and the ir~/ir + ratio 
provide complementary information about the symme- 
try energy. It has been argued that sub-threshold 77 pro- 
duction requires multiple nucleon-nucleon scattering over 
sufficiently long time [66[ , it may thus carry information 
mainly about the equilibrium phase of the reaction. At 
this stage, the following condition is satisfied between 
any two re gions of density (isospin asymmetry) pi(Si) 
andp 2 (5 2 ) [MEl 

E sym {pi)Si = E sym (p 2 )S2- (3) 

Then, the so-called isospin fractionation occurs, i.e., the 
density region with a higher symmetry energy E sym {p) 
will have a lower isospin asymmetry 5 and vice versa. 
Consequently, with a stiffer symmetry energy the isospin 
asymmetry 8 is lower in the high density regions where 
then more np (pn) collisions can occur to produce more 
77 mesons. While it is the opposite for the tt~ /tt + ratio. 

It is well known that the 7r~/7r + ratio is more sen- 
sitive to the symmetry energy at lower beam energies 
especially in the sub-threshold region where the mean- 
field dominates the reaction dynamics and has longer 
time to modify the momentum of nucleons. To create 
a 77 at sub-threshold energies requires even longer reac- 
tion time for both multiple collisions and the mean-field 
to act coherently in order to accumulate enough energy 
on the two colliding baryons or their resonances. Thus, 
it is interesting to examine the multiplicity of 77 pro- 
duction as a function of incident beam energy. Shown 
in Fig. [3] are the 77 multiplicity as a function of beam 
energy in inclusive Au+Au reactions with both x = 
and x = — 2. It is seen that the r\ multiplicity decreases 



rapidly with decreasing beam energy, especially for the 
soft symmetry energy. The 77 multiplicity saturates at 
an incident energy of about 10 GeV/nucleon. It is inter- 
esting to see that indeed the effect of nuclear symmetry 
energy is stronger in the deeper sub-threshold region as 
one expects. In fact, it is more sensitive to the sym- 
metry e nerg y than the tt~ /ir + ratio in the same energy 
region 109] although the pion yields are relatively more 
abundant. More detailed information about 77 produc- 
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FIG. 4: (Color online) Transverse momentum distributions of 
77 production in inclusive Au+Au reaction at incident beam 
energy of 0.2 GeV/nucleon with different nuclear symmetry 
energies and compression coefficients of symmetric matter. 

tion and its dependence on the symmetry energy can be 
obtained from studying its transverse momentum distri- 
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bution. Shown in Fig. [4] is such a distribution in terms 
of dMrj/dpt = dM v(pt-^Pt+dp t ) for Au+Au reactions at 

a beam energy of 0.2 GeV/nucleon with different values 
for the parameter x and the incompressibility K. It is 
seen that at a transverse momentum of about 200 ~ 300 
MeV/c, the effect of nuclear symmetry energy reaches 
its maxima of about a factor of 5-10. Changing the in- 
compressibility of nuclear matter from K=200 MeV to 
K=380 MeV seems to have little effect on the rj yields and 
its spectrum. While the larger (smaller) incompressibil- 
ity of nuclear matter causes smaller (larger) compression 
during the reaction, both the production and reabsorp- 
tion of rj are approximately equally affected. Moreover, 
as mentioned earlier, the p-n relative momentum is es- 
sentially determined by the gradient of the isovector po- 
tential with little influence from the isoscalar potential. 
Thus, the near-threshold rj production is sensitive to the 
symmetry energy but not the incompressibility of sym- 
metric nuclear matter. 



IV. PRODUCTION OF DARK U-BOSON FROM 
THE RARE DECAY OF v MESONS IN 
HEAVY-ION COLLISIONS 



the symmetry energy with x = can better reproduce 
the rj production data from the TAPS collaboration, we 
present here our estimates using the soft symmetry en- 
ergy with x = 0. Using the stiff symmetry energy with 
x=-2 would increase the maximal total U-boson produc- 
tion cross section by the same ratio shown in Fig. [3] 
Shown in Fig. [5] is the maximal total U-boson produc- 
tion cross section in Au+Au reaction at a beam energy 
from 0.2 to 10 GeV/nucleon. It is seen that the max- 
imum U-boson production cross section is in the 1-70 
jib range which is about 6-8 orders of magnitude higher 
than that estimated for the 7 + p reaction at the Jlab 
Eta Factory. However, similar to the situations at other 
facilities, it will be a big challenge to find in the dilep- 
tion spectrum in heavy-ion collisions the signature of the 
dark U-boson decay U —> £ + £~ because of the various 
backgrounds. Moreover, the predicted U-boson width is 
only about 10 eV which is much smaller than the resolu- 
tion of most detectors. Nevertheless, our estimate of the 
excitation function of U-boson production in heavy-ion 
collisions provide a useful reference for future U-boson 
search. 



V. CONCLUSIONS 
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In summary, within the transport model ART1.0 we 
studied the excitation function of rj yield in heavy-ion 
collisions both as a potential probe of the symmetry en- 
ergy and as a possible source of the dark U-boson. The 
yield of rj meson is found to be very sensitive to the 
density dependence of nuclear symmetry energy at sub- 
threshold energies. Using the branching ratio of the rare 
•q decay (rj — > jU) available in the literature, we esti- 
mated the maximum cross section for the U-boson pro- 
duction in heavy-ion collisions at a beam energy of 0.2 to 
10 GeV/nucleon, providing a useful reference for future 
U-boson search experimentally. 
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We now turn to estimating the maximum U-boson 
production cross section. For this purpose, we use the 
maximal U-boson-SM particle coupling constant e 2 = 
4. x 10~ 6 and BR(X -> Y + 7) ~ 0.7 |1 HH). Since 
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